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ANALYTIC STUDIES OF SOUND PRESSURES INSIDE THE DUCT 
OF DUCTED PROPELLERS 


By Joseph A. Drischler 
Langley Research Center 

SUIVIMARY 

The soimd-pressure field of a rotating ducted propeller in forward flight is analyzed 
by replacing the duct by an infinite coaxial circular cylinder and assuming that the blade- 
loading distribution associated with thrust and torque can be represented by a distribution 
of acoustic pressure doublets acting at the propeller disk. Trend studies are made to 
ascertain the effect on the pressure distribution inside the duct of variations in propeller 
loading, tip clearance, free-stream Mach number, tip Mach number, and hub-tip ratio. 

The effect on the pressure distribution of concentrating the loading at various radial posi- 
tions on the propeller blade is also investigated. 

One of the significant findings of the investigation is that the loading distribution 
(with the same total thrust and torque), the blade-tip clearance, and the free-stream 
Mach munber produce considerable variations in the radial sound-pressure distribution 
in the field near the propeller disk but produce insignificant variations at distances 
greater than about one duct radius from the propeller disk. Increasing the blade-tip 
clearance by only a few percent of the propeller diameter decreases the overall sound - 
pressure level by about 6 dB in the vicinity of the propeller tip. Increasing the free- 
stream Mach number generally increases the sound -pressure levels inside the duct, with 
the greatest increases occurring near the propeller disk and along the duct waU. At a 
Mach number of M = 0.7 the sound -pres sure levels aloi^ the duct wall are from 3 dB 
to 6 dB higher than for M = 0. 


INTRODUCTION 

At present, various types of aircraft being developed utilize ducted fans or pro- 
pellers. The presence of a duct increases the propulsion performance at low speeds and 
permits a reduction of propeller diameter. The presence of the duct presents a struc- 
tural problem not only for fans and propellers but also for compressors, and a knowledge 
of the sound field inside the duct is required to develop criteria for light-weight struc- 
tures of sufficient acoustic -fat%ue life, particularly in the critical region near the pro- 
peller disk. 
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Li reference 1 the sound pressures inside the duct of a ducted propeller were calcu- 
lated. The duct consisted of two infinite concentric circular cylinders, the inner cylinder 
representing the center body, and the wave equation was solved by the method of separa- 
tion of variables. The thrust and the drag (torque) forces on the propeller were combined 
to give a resultant force in the axial direction; no free-stream velocity through the duct 
was considered. In reference 2 the wave equation was also solved by the same technique. 
A free-stream velocity was introduced and the loading on the propeller was left arbitrary. 
In reference 3 Tyler and Sofrin also solved the wave equation by separation of variables 
and used assumed duct-velocity distributions to calculate the far -field sound pressures. 
However, the solutions were not combined to represent loadiigs on the ducted propeller 
nor was the pressure distribution in the duct derived. Neither reference 2 nor refer- 
ence 3 presents any numerical results on ducted propellers. The purpose of the present 
paper is to present a solution of the wave equation based on Fourier transform techniques 
and to present numerical results for the effects on the pressure distribution in the duct 
of variations in propeller loadii^, tip clearance, free-stream Mach number, tip A^ch 
number, and hub-tip ratio. The effect on pressure distribution of concentratii^ the 
loading at various radial positions on the propeller blade was also investigated. 

SYMBOLS 


a 

A 

■^nm ’ ®nm (^) 
b 
B 
c 

Hnm(r) 

^nm 

Jn.Yn 


k 


hub radius 
source strength 

undetermined functions needed to satisfy boundary conditions 
outer duct radius 
number of blades 
speed of sound 

quantity defined by equation (A23) 
quantity defined by equation (A12) 

Bessel functions of first kind and second kind, respectively 
nondimensional frequency parameter, ^ 




2 


M 


free-stream Mach number, V/c 


Mt 


n = n’B 


effective tip Mach number 
rotational tip Mach number 


n’ 

n’,n,m,p,q 

Q 

Q(p)/p 

Qnm 

ro,0,Zo 

r,e,z 

^cp 

Re 

SPL 

SPLqv 

t 

T 

V 

®n 


harmonic number 
summation indices 

pressure in n'th harmonic due to thrust and to torque, respectively 
total torque 

torque per imit nondimensional length 
quantity defined by equation (A9) 
cylindrical coordinates 

cylindrical coordinates normalized to outer duct radius b 

nondimensional center of pressure 

effective distance on span of propeller blade 

sound-pressure level, decibels 

overall sound-pressure level, decibels 

time 

total thrust 

free-stream velocity 

quantity defined by equation (A23) 


variable of integration 


Y 

6 

6(z) 

6 


^^nm 


Pq >^ ’^0 


P,0’,? 


P T(p) 

(7 

x(p) 


w 

S2 


notation defined by equation (A15) 

included angle of propeller blade considered in reference 1 

nondimensional propeller tip radius 

Dirac delta function 

imaginary part of k 

see equation (A6) 

mth root of = 0 

coordinates of source point 

nondimensional coordinates of source point 

thrust distribution, per unit nondimensional length 

hub- tip ratio, a/b 

velocity potential 

function in equation (1) denoting loading distribution 
<»■ 

PQ 

frequency of n'th harmonic 
propeller angular velocity 

ANALYSIS 


The sound-pressure field of a rotating ducted propeller is analyzed on the assump- 
tion that the pressure distribution over the propeller, associated with the thrust and the 
torque, can be represented by a distribution of rotating acoustic pressure doublets acting 
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at the propeller disk and subject to uniform rectilinear motion. This work is a general- 
ization of the work done by Garrick and Watkins, for a free propeller reported in refer- 
ence 4. It is further assumed that the pressure field near the propeller can be adequately 
predicted by analyzing the field between two infinite coaxial circular cylinders, the iimer 
cylinder representing the center body, usually small in ducted propellers and large in 
prop-fans or compressors. The analysis does not take into account reflections, such as 
would occur at inlets; however, these reflections are generally small for ducted fans, 
except possibly very near the inlet, because of the cutoff phenomena for subsonic spinning 
modes. 

Figure 1 is a schematic of the ducted-propeller disk and the coordinate system 
associated with the mathematical model used in the present analysis. The entire analysis 
is presented in the appendix where is derived the velocity potential or Green's function 
which satisfies the wave equation for a pulsating source rotating with angular velocity SI 
and moving in the negative z-direction with velocity V. By properly representing the 
forces due to thrust and torque (or drag) these sources are superimposed to yield the 
pressure distribution inside the duct due to B blades rotating with ai^ular velocity SI. 

RESULTS AND DISCUSSION 

Four distributions of thrust and torque are utilized to determine their effect on the 
pressure distribution inside the duct. The distributions are chosen so as to have the 
center of pressure between the 0.6- and 0.8-spanwise stations and to be representative of 
the loadings of most ducted propellers. All distributions are normalized so as to yield 
the same total thrust T and total torque Q; that is, 

T = B p T(p) dp 


and 




where 6 and a are 
By letting x(p) 


the nondimensional tip and 

denote either or 

T 


hub radius, respectively. 

^ the four cases considered are 

pQ 
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In all cases 6 is chosen to equal 1.0 except for some studies of tip clearance in 
Case IV which are discussed subsequently. Case I assumes that all the thrust and all the 
torque are concentrated at an effective radius equal to 0.8 of the duct radius. This same 
assumption is discussed in reference 4 for a free propeller. Case II represents a para- 
bolic loading distribution and Case HI is loaded according to J8(^lP) which has its 
center of pressure at approximately 0.75 the blade radius with ^ 8 (^ 1 ) = ^ where is 
the first root of J8(x) = 0. Case IV represents a linear loading distribution. For 
comparison, the four basic types of propeller loading are presented in figure 2 for a hub- 
tip ratio of 0.2. The values for the centers of pressure are also given. 

Values of thrust, torque, blade radius and number, rotational speed, and blade angle 
for a particular ducted propeller were obtained from reference 1 and are given with other 
pertinent data in table I. These data are used throughout the present report for all 
numerical studies lonless otherwise stated. Measurements and calculations were made in 
the U.S. Customary units. The values are presented herein in the International System of 
Units (SI) with the equivalent values given parenthetically in the U.S. Customary Units. 

Radial Pressure Distribution 

In figures 3 to 6 the overall root-mean- square (rms) radial pressure distribution 
is presented for M = 0 at various values of z, the distance from the propeller disk. 

The nondimensional radial coordinate varies from r = 0.2 at the hub to r = 1.0 at the 


6 



duct wall. It can be seen in figure 3 that large differences occur in overall sound- 
pressure level SPLoy among all four cases, the greatest difference being in Case I 
at r = 0.8 where the loading is concentrated. Farther downstream (figs. 4 to 6) the 
effect of the loading distribution becomes less apparent and, although not presented, at a 
station approximately one duct radius from the propeller disk, differences in the SPLqv 
among the four cases become indistii^uishable. This is as expected since examination of 
equations (A24) and (A25) indicates that the pressure drops off exponentially as the dis- 
tance from the propeller disk increases. At z = 0 the exponent of the exponential 
vanishes and more terms are needed to define the SPLpy. A further indication of this 
trend may be noted in figure 7 where the SPLqv along the duct wall (r = 1.0) is plotted 
against z. At z = 1.0, although not shown, further calculations indicate that the first 
term dominates and the values of SPLoy in all four cases converge to within 1 dB. 

Whereas figures 3 to 7 show the SPLoy, figures 8 to 10, in which sovind-pressure 
level SPL is plotted against harmonic number n' , show the contribution of each 
harmonic to the values presented in figure 7 for SPLqv along the duct wall at r = 1.0. 

In figure 8, which represents the SPL in the propeller disk, it appears that all 
harmonics shown contribute significantly to the SPLoy, ^ indication that the 20 terms 
used in the series and as mentioned in the appendix were insufficient to determine the 
pressure in the plane at z = 0. This is due to the fact the exponent of the exponential is 
zero in the propeller disk whereas, downstream of the disk (at z = 0.06 and at z = 0.1 
in figs. 9 and 10, respectively), since the series converges rapidly, the contribution from 
the higher harmonics falls off with increasing harmonic number. It may also be noted in 
figures 9 and 10 that downstream of the propeller disk, in the first two harmonics, values 
of SPL agree to within 3 dB for all four cases; however, in the higher harmonics there 
are large differences which are believed to be due to the details of the different loading 
distributions. 


Effect of Tip Clearance 

In reference 1 are presented the results of some experimental measurements of the 
effect of tip clearance on the overall sound-pressure level in the duct. These measure- 
ments are taken on the propeller for which the parameters are given in table I except 
that values of thrust and torque are replaced by T = 5224.6 N (1174.6 Ibf) and 
Q = 1783 J (1315 ft-lbf) to correspond with the experimental results presented in 
figure 17 of reference 1. 

In the present paper only Case IV (linear loading) is considered to illustrate the 
effect of tip clearance. The clearances used herein are the same as those of reference 1, 
namely, 0.95 cm (3/8 in.) and 3.81 cm (1.5 in.), which correspond to values of nondimen- 
sional blade tip radius equal to 0.991 and 0.964, respectively. The effect on the SPLqv 
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across the propeller disk is shown in figure 11 and along the duct wall (r = 1.0), in fig- 
ure 12. Here it can be seen (fig. 11) that the tip clearance is effective in reducing 
the SPLov over only the outer 5 percent of the disk and has a negligible effect on the 
remaining portion of the disk. However, as mentioned previously, and in the appendix, the 
results for z = 0 (at the propeller disk) are questionable. Along the duct wall (fig. 12) 
the effect seems to extend from z = 0 to approximately z = 0.2. This reduction in 
pressure could be significant in the structviral design of ducted propellers; however, the 
loss of performance may be more important in determining the design criteria. 

Also plotted in figure 12 are the experimental data obtained from figure 17 of refer- 
ence 1. The symbols with ticks represent data points taken upstream of the propeller 
disk and are included in figure 12 since the present analysis indicates equal attenuation 
in both directions. As can be seen in figure 12, the measured pressures upstream of the 
disk are greater than those downstream and the differences due to tip clearance are much 
greater than those indicated by the theory. The overall agreement between experimental 
and theoretical results is felt to be poor and no explanation of these differences can be 
given at the present time. 

Effects of Free-Stream Mach Nimiber 

Two different free-stream Mach numbers were considered for comparison with the 
static case M = 0, namely, M = 0.3 and M = 0.7, which yield effective Mach numbers 
of Me = 0.73 and Mg = 0.97, respectively. Only Case H (parabolic loading) was consid- 
ered and the overall sound-pressure level is presented for M = 0 and M = 0.7 in fig- 
ures 13 and 14. The restilts for M = 0.3 are not shown since they are essentially the 
same as those for M = 0. In figure 13, the SPLqv at the duct wall shows an increase 
of between 3 dB and 6 dB due to Mach number effects as the station moves downstream 
(or upstream) from the propeller disk at z = 0 to z = 0.4. In figure 14 it can be seen 
that the difference decreases as the station aproaches the hub; for z = 0.1 the effect of 
Mach number is negligible. It may be noted here, and is dis cussed more fully in the 
appendix that, for effective Mach numbers Mg = ^M^ + Mj2 < 1.0, Mj being the rota- 
tional tip Mach number, the pressure attenuates equally in both upstream and downstream 
directions. 


Effect of Rotational Tip Mach Number 

For the case of linear spanwise loading (Case IV), the effect of rotational tip Mach 
number is investigated for two free-stream conditions, M = 0 and M = 0.7. These 
results are presented in figures 15 and 16, which show the overall sound-pressure level 
along the duct wall for several rotational tip Mach nu mbers betw een 0.6 to 2.0. Here it 
can be seen that for an effective Mach number Mq = y^M2 + M^2 < i.o the pressure 
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decreases exponentially, an Indication of attenuation of the pressure waves, whereas 
for Me >1.0 the waves propagate in both directions with nearly equal amplitude all 
along the duct wall. The criterion for propagation, as mentioned previously, is discussed 
in the appendix; however, it may be illustrated by examining the exponential in equa- 
tions (A24) and (A25). Wlien 1® ^ quantity; that is, when 

> (l - 


or 


(nUtf > (l - M2)Mnm^ 

the contribution from that term in the series causes propagation; otherwise the wave 
decays exponentially. For example, for M = 0 and Mj. = 1.2, Pnm luust be less 
than 1.2n, where n = n'B, in order for a wave to propagate in the duct. The tables of 
reference 5 show this inequality to be satisfied when n § 9 (n' =3). This is further 
illustrated in figure 17 which gives the harmonic content of duct-wall pressure 
for M = 0 and Mj = 1.2. It may be noted that the pressure harmonics attenuate for 
n' = 1 and n' = 2 whereas the higher harmonics propagate. 

Results are shown in figures 15 to 17 to indicate how various harmonics can decay 
or propagate depending on the effective Mach number. Although effective Mach numbers 
greater than one are physically possible for configurations such as prop-fans, the Mach 
numbers presented herein are higher than the operating range of the particular configura- 
tion considered in reference 1. Also, since figure 17 indicates that the propagating waves 
produce sound-pressure levels of nearly the same magnitude (at least to the tenth 
harmonic), it would appear inappropriate to sum algebraically the squares of the pressures 
in each mode. A further investigation shoxdd be undertaken to determine the manner in 
which the individual modes should be combined. 

Effect of Hub-Tip Ratio 

The effect on the overall sound-pressure level of changing the hub-tip ratio 
from cr = 0.2 to a = 0.5 is illustrated in figures 18 to 20 for the linearly loaded blade 
(Case IV). In figure 18 the radial variation in SPLqv is presented for stations z = 0.2 
and z = 0.4. Here it can be seen that concentrating the loadii^ over the outer half of the 
blade resvilts in an increase in SPLqv along the duct wall of approximately 2 dB 
to 3 dB over that for the more lightly loaded blade even though the rotation speed, 
thrust, and torque are the same for both cases. This same trend is also evident in fig- 
ure 19, where the SPLqv is given along the duct wall and the difference in decibels 
remains essentially constant as the station moves downstream. Also, in figure 20 where 
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the harmonic content of the overall SPL at z = 0.1 and r = 1.0 is shown, a differ- 
ence of approximately 2 dB exists for all 10 harmonics. 

Effect of Varying Effective Radius of Concentrated Load 

The concept of concentrating the thrust and the torque at some effective distance Rg 
on the span of the blade was used by Garrick and Watkins in reference 4. The suggested 
value of Rg = 0.8 (ref. 4) is used herein, together with values of Rg = 0.6 
and Rg = 0.7. The results of these calculations are presented in figures 21 to 23; 
figure 21 shows the variation in overall sound-pressure level in the radial direction and 
figure 22, the variation along the duct wall. Figure 23 shows the harmonic content of the 
overall sound pressure along the duct wall at z = 0.1 and r = 1.0. It can be seen from 
figures 21 to 23 that the effective-radius concept leads to large differences in sound- 
pressure level near the propeller disk. By comparing figures 7 and 21 it can be seen 
that, even though the concentrated-loading theory yields pressure distributions which are 
essentially independent of the loading as the station moves downstream in the duct, it 
should be used with caution near the propeller disk. 

Comparison of Present Results With Reference 1 

In reference 1, the wave equation for a propeller in an infinite coaxial circular 
cylinder is solved by separation of variables. The pressure on the propeller blade is 
assumed to be the resultant of two forces, the thrust acting in the axial direction and the 
drag acting in the circumferential direction at the center of pressure rgp. The resultant 
of these two vectors, that is. 



was assumed to be distributed uniformly over the blade (a sector of angle y) and to act in 
the axial direction. The absolute value of the Fourier coefficient of the n*th harmonic for 
the loading can be represented by 


fv \ 

Po 

slnSllS 


irb2(l - a^) 

n'By 

2 


This loading differs from the loadir^ considered in the present paper where the thrust and 
the torque are considered concentrated aloi^ radial lines and thus yielding a comparable 
equation 



( 3 ) 
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r 


where Pq* is representative of either the thrust force T or the drag force Q/r^p. 
Equations (2) and (3) have the same form as y — 0. 

Circumferential integration over each blade considered in reference 1 results in a 
linear pressure variation in the radial direction and thus is comparable to case IV 

x(p) = — considered herein. By using the notation of the present paper, the n'th 
1 - cr^J 

larmonic of the pressure distribution presented in reference 1 can be written in the form 


(Prms)n,- 


8\^Pq sin^ 
nyb^(l - (T^) 


, 1/2 


2 Hnm(r) ofn J p ^jp 


m=l' 


( 4 ) 


and 


Prms 

where, as before, n = n'B. A comparison of equations (A24) and (4) shows that with the 
exception of the loading, the pressure distribution in the duct is given, in the present 
analysis, as the sum of the radial modes whereas, in reference 1 , the pressure distribu- 
tion is the square root of the sum of the squares of the radial modes. 



Figure 24 shows the harmonic content of the overall SPL for some calculations 
based on equations (4) and (A24), together with experimental data from figure 31 of 
reference 1. The data were obtained from reference 1 and are for z = 0.104, r = 1.0, 
and M = 0. Other pertinent quantities are presented in table I, also obtained from 
reference 1. In figure 24, it can be seen that the results from both theories and the 
experimental results are in fair agreement at least to the sixth harmonic. The differ- 
ences thereafter are believed to be due to the differences in the details of the loading. 

As noted in equation (2) the contribution from the n*th harmonic of the pressure loading 
is zero when ^ = 77 or when n' = Utilizing the constants y = 0.27 rad 

and B = 3 from table I yields a value of n* = 7.75. Examination of figure 24 indicates 
a great reduction in SPL for a value of n’ = 8 which, as mentioned previously, may be 
due to the differences in the details of the loading used in the present paper and that of 
reference 1 . 


SUMMARY OF RESULTS 

The sound -pressure level inside the duct of a ducted propeller is analyzed by 
replacing the duct by an infinite coaxial circular cylinder and assuming that the blade- 
loading distribution associated with thrust and torque can be represented by a distribution 
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of acoustic pressure doublets acting at the propeller disk. Trend studies are made to 
determine the effect on the pressure distribution inside the duct of variations in propeller 
loading, tip clearance, free-stream Mach number, tip Mach number, and hub-tip ratio. 

The effect on pressure distribution of concentrating the loading at various radial positions 
on the propeller blade is also investigated. 

The following results are observed: 

1. The radial pressure distribution at least one radius downstream of the propeller 
disk is essentially independent of the blade loading; however, there is a significant varia- 
tion in the field near the propeller disk. 

2. Increasing the blade-tip clearance by only a few percent of the propeller diameter 
has a significant effect in reducing the overall sound-pressure level along the duct wall by 
about 6 dB in the field near the propeller disk; however, at approximately one radius from 
the disk, the effect of tip clearance is negligible. 

3. For the angular velocity considered, free-stream Mach numbers equal to or less 
than 0.3 have very little effect on the pressure distribution whereas, for M = 0.7, the 
sound-pressure levels along the duct wall are from 3 dB to 6 dB higher than for M = 0. 

4. For effective tip Mach numbers Mg < 1 the sound waves attenuate in the duct 
whereas for Me > 1 the waves propagate, the sound-pressure level along the duct wall 
remaining essentially constant and changing only near the propeller disk. 

5. For the case with linear distribution and the same total load, increasing the hub- 
tip ratio from a = 0.2 to 0 = 0.5 results in an increase in sound-pressure level of 
approximately 2 dB to 3 dB throughout the duct. 

6. Changing the location of the effective radius Re of the concentrated load results 
in significantly different sound-pressure levels near the propeller disk. As Rg is 
increased the pressures on the outer wall increase and those on the inner wall decrease; 
however, as pointed out previously in item (1), as the distance downstream increases, the 
pressure distribution across the duct becomes independent of the load distribution. 

7. As the effective radius is decreased, the sound-pressure levels in the higher 
harmonics rapidly decrease, so that it is necessary to calculate only a few harmonics to 
determine the overall sound -pressure level. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Virginia, June 23, 1971. 
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APPENDIX 


SOUND-PRESSURE LEVEL INSIDE AN INFINITE ANNULAR DUCT 

The purpose of the appendix is to derive the pressures due to thrust and torque 
in an annular duct containing uniform subsonic or supersonic flow. First, expressions 
are derived to give the pressures due to unit pulsating forces acting in the axial direction 
to simulate the thrust and in the tangential direction to simulate the torque, or drag. By 
superposition, the forces are combined in such a way as to represent the thrust and the 
torque distribution on a propeller in the duct. A Fourier time expansion of the thrust 
and the torque is made to account for the rotating blades. This procedure allows each 
blade harmonic to be treated as a distribution of pulsating forces. 


Derivation of Pressure Due to Concentrated Pulsating Force 
in Duct of Infinite Length 

The velocity potential due to a pulsating source of strength A in an annular duct 
with free-stream velocity V directed along the positive Zg-axis is governed by the 
wave equation 

’^*0 = 3 *0* (■^0 - Po) 5 ("0 - «o) - »’) Wi) 

with boundary conditions such that the velocity normal to the, duct walls be zero; that is. 



where rg = a, the inner or hub radius, and rg = b, the outer or wall radius. In equa- 
tion (Al), rg, 0, and Zg are the coordinates of the field point and pg, 0’, and ?g 
those of the source point. Figure 1(a) is a schematic of the propeller disk and the coor- 
dinate system associated with equation (Al). Figure 1(b) is the mathematical model used 
in the present analysis. The coordinates are normalized with respect to the outer-duct 
radius where r = cr (o' = and r = 1.0 correspond to the hub and the outer wall, 
respectively, and the duct is of infinite length (-oo < z < °o). 

By letting <^g = <pe and normalizing all coordinates with respect to b, equa- 
tion (Al) becomes 
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APPENDIX - Continued 


and 



where M = ^ and k = In the following analysis it is assumed that k is a complex 
quantity with a positive imaginary part. This assumption is equivalent to introducing a 

damping term — ^ into the wave equation (Al). As is seen later, in connection with 
equation (A12), this device moves the singularities off the path of integration and simpli- 
fies the evaluation of certain integrals. The final solution is obtained by lettii^ the 
imaginary part of k approach zero. A further discussion of this can be found, for 
example, in reference 6. 

The solution of equation (A2) is assumed to have the form 


CO 

" Z ^n(/3) + Yn(/ir) Bn(^)J d/3 

n=-oo 


(A4) 


where Jjj(/u,r) and Yjj(/xr) are Bessel functions of the first and second kind, respec- 
tively, and An(j3) and are arbitrary functions of /3 which must be chosen so as 

to satisfy the boundary conditions of equation (A3) . The functions A^O) and Bjj(^) 
are determined by substituting the expression in the square bracket in equation (A4) into 
equation (A3). Thus, equation (A4) becomes 


n=-oo m=l V ' 

in which the eigenfunctions given by 

?n(^nm^) = Jn(^^nm’^) “ '^n’(^nm) Yj^^Mnm^) 


(A5) 


(A6) 


and iJ-nm is the mth eigenvalue of the characteristic equation 


?n'(Mnm^)=0 (A7) 

where the primes are used in equations (A5), (A6), and (A7) to indicate differentiation of 
the function with respect to its argument. Reference 5 contains tables listing values of 
ju, ^^ for n = 0 to n = 25 and m = 1 to m = 10 with values of a ranging from 0 
to 0.9 in steps of 0.1. It may also be noted that 
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APPENDIX - Continued 


E 



'^n' (Mnm) 


where E is the notation used by Tyler and S'ofrin in reference 3. 

Substitution of equation (A5) into equation (A2) yields 

n=-oo m=l 

= - ^ 6(r - p) 5(z - C) 6(0 - 0’) (A8) 

Using the orthogonal relations 


r2^gi(n-n‘)0 de = <f° 

‘^0 |^ 27 T 


(n * n’) 
(n = n') 


and 


> (A9) 


^nKm^) ?n(f^nm'r) r dr = 0 


Qnm “ ^ 


M'nm 


^Un^(Mnmr) 


r=l 


r=cr 


* m') 
(m = m') 


ogether with the inverse Fourier transform with respect to /3 yields the expression 


■"^n' (^nm) ^ ?n (/^nmP) 


*nm 


4’r2Qnmp(l - m 2) - 2/3kM + 

K equation (AlO) is substituted into equation (A5), the result can be written in the form 


(AlO) 


<p(r,p,e,z) = - 2 ^ 


00 - Agin0 


n=-oo m=l 


^n(^nm^) ^n(^nmP) ^ 

4ir2bQ, 


nm 


(All) 


«nm 
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APPENDIX - Continued 


where 


*nm 


=i: 




(l - m2)/s 2 . 2kM/3 + Mnm^ - 


(A12) 


and 


e = e - e’ 


z = z - ? 

Here <f>(r,p,9,z) is the Green's function which represents the velocity potential at r,6,z 
due to a concentrated source of strength Ae~^^^ located at p,6' in a fluid with free- 
stream velocity V. 


Evaluation of Inm> Equation (A12) 

The integral I nm can be evaluated by the theory of residues. To accomplish this 
it is necessary to investigate the zeros of the denominator of the integrand, namely 




M(k + ie) ± - (l - M^)|Linm^ + 2ike 


1 - 


where e is the positive imaginary part of k and, as mentioned previously, has the 
effect of moving the poles off the path of integration. For M > 1.0 it can be shown that 
both poles lie in the lower half of the ^-plane. In the limit as 6 — 0 the path of inte- 
gration passes above both poles which are located at 


(A13) 



and the value of the integral is, for M >1.0, 



M>1.0 


_2^gikMz/(M2-l) z/k2 +(m2 - 
, - -- -= Sin-1 

/k2 + (M2-l)Mnm^ ^ 


(z > 0) 




(A14) 


= 0 


(z < 0^ 
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APPENDIX - Continued 


For M < 1.0 one pole lies in the upper half of the 0-plane and the other in the lower half 
of the 0-plane. In the limit as e — 0 the path of integration passes between the poles 
which are located at 




(A15) 


where the upper sign applies when the pole is in the upper half of the plane and the lower 
si gn when the pole is i n the lower half of the plane. By adopting the notation 

+ y - (l - (l - equation (A15) can be rewritten 




(A16) 


where the sign convention defined in equation (A15) applies. With use of this definition 
there can easily be obtained, for M < 1.0 _ 


Wm<1.0 


- (i - 


(z<0) 


/k2 - (l - 


(z>0) 


Substituting either equation (A14) for M > 1.0 or equation (A17) for M < 1.0 into 
equation (All) yields the velocity potential for a concentrated source of strength Ae“^^^ 
By the proper choice of A it is shown in reference 4 that, for a force 

T^Po>®'»?o) Pq ‘^Po 

acting in the z-direction, the pressure at any field point Tq,6,Zq can be obtained by 
differentiating the velocity potential with respect to Zq and replacing A 
by T(pQ,e',?Q) Pq dp(j d0' so that 

n_ = T/n-.fl'.r.N n- rtn- rffl* 1 . 1 ^ 




or in nondimensional quantities 


p.p = bT(p,e’,?)pdpde'i-^ 


(A19) 


I 


APPENDIX - Continued 

where T(p,0’,O is the thrust per unit area. Similarly, if A is related to 


Q(Po'^'go) 
fo 


Pq ^Pq 

a force in the tangential direction, the pressure can be expressed as 


PQ = Q(p.P%?)<sp<ie’^|^ (A20) 

and q(p, 0',?) is the torque per unit area. The problem now remains to express (in a 
Fourier series) the thrust and the torque as a distribation of concentrated forces of the 
form Ae”^^*^ so that each term of the expansion is of the form for which equation (All) 
was developed. 


Fourier Expansion of Thrust and Torque Distribution on Rotating Blades 

Each element of the propeller blade is acted on by the surface-pressure distribution 
which may be resolved into a thrust force in the direction of the Z-axis and a force asso- 
ciated with the torque acting in the circumferential direction. A description of these 
forces, including a discussion of the effect of blade thickness on the sound pressure, and 
their Fourier expansions are given in reference 4. The present analysis to this point has 
been general in regard to the distribution of thrust ana torque; now the analysis is 
specialized to the case of forces on a propeller which lies in the plane where z = 0 and 
rotates with angular velocity Ji. Further, the forces are assumed to be distributed 
radially along B blades ^the angle between the blades being and to be concentrated 
circumferentially at 


0’ = m + 


2irp 


(p = 0,1,2, . . . (B-1)) 


This can be represented by 


B-1 


p=0 ^ 


(A21) 




where \j/ is representative of either the thrust T^Pq, 0’^ or the drag force Q( 
per unit area and represents the thrust or the drag force per unit lei^th. This 

series when expanded in its Fourier form becomes 


f' = M»o) I I ^ ^ 


(A22) 


p=0 q=-o 
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APPENDIX - Continued 
n 

from which it can be seen that k = - and each term in the series is of the form 

f ^ 1 . c 

Ae” for which equation (All) was developed. It now remains to substitute the appro- 
priate expressions in equations (A19) and (A20) for prp and Pq, respectively, and to 
integrate over the propeller disk to obtain the total pressure due to thrust and torque. 


Pressure Due to Thrust and Torque 

Since supersonic flow through a ducted propeller is of only academic interest, 
equations that apply to subsonic flow are considered from here on; however, supersonic tip 
Mach numbers are not excluded. Substituting the appropriate ejq>ression for (p from 
equation (All) and for Inm from equation (A17) into equation (A19) and integrating over 
the propeller disk gives the pressure p>j> due to the thrust distribution T(p) 


'T- 2, III ^ [ i, ’ ^ 

n=-oo m=l p=0 


where 


a 


n 


Mk ±- 






2 

nm 


1 - m2 


Hnm(r) 



-‘«n( 


1 



(i-m2) 

1 2 

' f^nm 



> (A23) 




Qnm is defined by eqioation (A9), and a and 6 are nondimensional distances to the 
hub and the tip, respectively. By using the equalities 


g-i(n-q) 0 ’ ^ 0 . = [ ^ 


and 


B-1 . 2irp r 

y 

^0 V 


(o 


=„) 

*n) 


(n = n' b) 
(n n’ b) 


19 


APPENMX - Continued 


choosing the reference angle 6 = 0 (since all other values of 9 have, with the excep- 
tion of phasing, the same time history), dropping the time factor, and combining the 
±n terms, it can be shown that the n'th harmonic of the pressure due to thrust is 

op g 

(Pt)^, = ^ ^ o^e-i«nZH^(r)£ Bp T(p) dp (A24) 

Similarly, the n'th harmonic of the pressure due to torque is 

(pQ)n’ ~ X ®^nmW J ^n(^nmP) p^ dp (A25) 

“ m=l ^ 

where = 1 for n' = 0 and = 2 for n' 0 and the root -mean-square (rms) 
pressure in the n'th harmonic is 


(Prms)„. = 


±_ 

1/2 


N„. * H. 


The overall rms pressure is given by 


Prms 



and the overall soimd-pressure level in decibels is given by 


(A26) 


(A27) 


SPLqv = 128 + 20 logjQ Prms (A28) 

if Prins is inlb/ft2. 

Equations (A24) to (A28) have been programed for the four thrust and torque distri- 
butions given in equation (1). The first 10 harmonics are usually presented for the condi- 
tions considered and the summation on m contains 20 terms. For values of z > 0.1, 

5 terms appear sufficient for convergence; however, for z = 0, althovigh the series 
appears to converge within 12 terms, it is questionable whether 20 terms are sufficient to 
give accurate results. 
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APPENDIX - Concluded 


Criteria for Attenuation or Propagation of Wave in Duct 


Examination of equations (A24) and (A25) shows that, when Ofn is complex, certain 
terms in the series indicate oscillations or propagation downstream whereas when afjj 
is imaginary, they indicate exponential decay as z increases. This leads to an investi- 
gation of whether - (l - is greater than or less than zero. By substituting 

nMj. for k and since the expression can be rewritten as 



(A29) 


where, when the upper inequality in equation (A29) is satisfied, propagation is indicated 
and, when the lower inequality is satisfied, decay or attenuation occurs. The quantities 
to be examined are Me 1.0 and M ^ 1.0 For Mg < 1.0, which implies that M < 1.0, 
the left-hand side of equation (A29) becomes negative and the lower inequality is satisfied 
since, as indicated in reference 7, f^nm waves attenuate. For Mg > 1.0, 

two alternatives exist, that is, M > 1.0 or M < 1.0. For M > 1.0, the right-hand side 
of equation (A29) becomes negative, the upper inequality is satisfied, and all waves 
propagate. For M <1.0, when the approximate asymptotic expansion of the mth zero 
/inm ® considered it can be seen that, for a fixed m, a value can 

.^always be assigned to n for which the upper inequality is satisfied and that all circum- 
ferential modes for which the mode number is equal to or greater than this value of n 
propagate in the duct. It may also be noted (ref. 7) that, for large values of n, 
the eigenvalues are approximately the zeros of Yn’ (j'^nm) ~ first eigenvalue 

being 


Mn,l “ n + 1 . 82 n^/^ 

When this expression is substituted into equation (A29) it can be seen that Mg - 1 is of 
the order of 1/4^/^, again an indication that for a large enough value of n the upper 
inequality is satisfied and all modes with mode numbers greater than or equal to this 
value of n also propagate. This is also illustrated in figure 17 and mentioned in the 
text, that for M^ = Mg = 1.2 and M = 0, the first two modes attenuate whereas the 
remaining modes propagate. It therefore can be concluded that for Mg > 1.0 an infinite 
number of modes, but not all, propagate and for Me <1.0 all modes attenuate. 
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TABLE I.- DATA FOR DUCTED PROPELLER 
TESTED IN REFERENCE 1 

T 6921 N (1556 Ibf) 

Q 2766 J (2040 ft-lbf) 

S2 2000 rpm = 209 rad/sec 

b 1.067 m (3.5 ft) 

B 3 

y 0.27 rad 

<T 0.2 

c 335 m/sec (1100 ft/sec) 
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Figure 7.- Variation of sound-pressure level along 
duct wall at r = 1.0. M = 0. 
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Figure 8.- Variation of sound-pressure level with harmonic 
number at station z = 0. r = 1.0; M = 0. 
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Figure 9.- Variation of sound-pressure level with harmonic 
number at station z = 0.06. r = 1.0; M = 0. 
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Figure 10.- Variation of sound-pressure level with harmonic 
number at station z = 0.1. r = 1.0; M = 0. 
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Figure 11.- Effect of tip clearance on sound-pressure level in 
radial direction at z = 0. M = 0; T = 5224.6 N (1174.6 Ibf); 
Q = 1783 J (1315 ft -Ibf); Case IV. 
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Figure 12.- Effect of tip clearance on sound-pressure level along 
duct wall at r = 1.0. M = 0; T = 5224.6 N (1174.6 Ibf); 

Q = 1783 J (1315 ft-lbf); Case IV. 
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Figure 13.- Effect of free-stream Mach number (M = 0 and M = 0.7) on 
soimd-pressure level in downstream direction for r = 1. Case II. 
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Figure 14.- Effect of free-stream Mach number (M = 0 and M = 0.7) on radial 
pressure distribution at two downstream locations, z = 0 and z = 0.1. Case II. 
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Figure 15.- Effect of rotational tip Mach number on sound-pressure level for M = 0. r = 1.0; Case IV. 






r 


Figure 18.- Radial variation of sound-pressure level for values of hub-tip 
ratio of o = 0.2 and a = 0.5 at z = 0.2 and z = 0.4. M = 0; Case IV. 
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Figure 19.- Variation of sound-pressure level at duct wall for values of 
hub-tip ratio of ct = 0.2 and a = 0.5. M = 0; r = 1.0: Case IV. 




20.- Variation of sound-pressure level with harmonic number for values 
of hub-tip ratio of cr = 0.2 and a = 0.5. z = 0.1; M = 0; r = 1.0; Case IV. 
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(b) z = 0.4. 

Figure 21.- Radial variation of sound-pressure level for values of effective 
radius Rg = 0.8, Rg = 0.7, and Rg = 0.6 at z = 0.2 and z = 0.4. 

M = 0; Case I. 



Figure 22.- Variation of sound-pressure level along duct wall for 
values of effective radius Rg = 0.8, Rg = 0.7, and Rg = 0.6. 
M = 0; r = 1.0; Case I. 
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Figure 23.- Variation of sound- pressure level with harmonic 
number for values of effective radius Rg = 0.8, Rg = 0.7, 
and Rg = 0.6. z = 0.1; M = 0; r = 1.0; Case I. 




Exp. (Ref. 1) 

Theory (Ref. 1) 

Present theory (Case IV) 


Figure 24.- Comparison of experimental and theoretical results shown in 
figure 31 of reference 1 with results of present analysis, z = 0.104; 
M = 0; r = 1.0; Case IV. 
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